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Abstract. The equations of motion of two point masses in harmonic coordinates are derived 
through the third post-Newtonian (3PN) approximation. The problem of self-field regularization 
(necessary for removing the divergent self-field of point particles) is dealt with in two separate steps. 
In a first step the extended Hadamard regularization is applied, resulting in equations of motion which 
are complete at the 3PN order, except for the occurence of one and only one unknown parameter. In 
a second step the dimensional regularization (in d dimensions) is used as a powerful argument for 
fixing the value of this parameter, thereby completing the 3-dimensional Hadamard-regularization 
result. The complete equations of motion and associated energy at the 3PN order are given in the 
case of circular orbits. 
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1. Introduction 

The third post-Newtonian approximation of general relativity (in short 3PN),^ became 
famous in recent years because of its frightening or, depending on one's state of mind, 
fascinating intricacy. In particular, the study of this approximation reached a somewhat 
paroxysmic stage when it was realized that the usual self-field regularization, based on 
Hadamard's concept of "partie finie" [2,3], although having proved to be very efficient 
up to the 2PN order, fails to provide a complete answer to the problem at the 3PN order. 
Indeed it seems to inexorably yield the appearance of some numerical coefficients which 
cannot be determined within the regularization. 

Working at such a high approximation level as the 3PN one does not represent a purely 
academic exercise. The current network of laser-interferometric gravitational-wave detec- 



*To appear in the Proceedings of the International Conference on Gravitation and Cosmology 
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^Following the standard practice [1], we refer to nPN as the terms of the order of in the 

equations of motion, relatively to the Newtonian acceleration. 
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tors (notably the large-scale ones: LIGO and VIRGO) will soon make possible the study of 
the inspiral and coalescence of binary systems of neutron stars and black holes. To extract 
useful information from the gravitational waves, theoretical general-relativistic waveforms 
are used as templates in these experiments, and it has been demonstrated that these must 
be extremely accurate, which means probably as accurate as the 3PN approximation [4-6]. 
To construct the 3PN templates one needs to control both the binary's equations of motion, 
at the 3PN order relatively to the Newtonian acceleration, and the gravitational radiation 
field, also consistent at 3PN order but with respect to the famous Einstein quadrupole for- 
mula, corresponding to the "Newtonian" order in the waveform. 

In this paper we focus our attention on the problem of motion of a point mass binary sys- 
tem. The undetermined parameters which appear, due to Hadamard self-field regulariza- 
tion, are several, but, in fact, once one has invoked physical arguments to compute some of 
them, it remains one only one unknown coefficient, the so-called "static" ambiguity param- 
eter ujs in the 3PN Hamiltonian in ADM coordinates [7,8], or, equivalently, the parameter 
denoted A in the 3PN equations of motion in harmonic coordinates [9,10]. [We mean by 
physical arguments the requirement of invariance under global Poincare transformations, 
and the demand that the equations of motion should be derivable from a Lagrangian (ne- 
glecting the 2.5PN radiation reaction term).] These parameters are related to each other by 
[9,11,12] 

^ = -11'"'^ - 3080- 

On the other hand, concerning the radiation field, three other parameters, ^, k and C, com- 
ing from the Hadamard regularization of the 3PN quadrupole moment, appear [13]. There 
is, however, a single parameter which enters the orbital phase of inspiralling compact bi- 
naries, in the form of a Hnear combination of ^ = ^ -|- 2k; -|- C and A. [Notice that A 
enters the radiation field because of time differentiations of the 3PN quadrupole moment 
and replacement of the accelerations by the 3PN equations of motion.] 

The regularization "ambiguities", say or A, are not real physical ambiguities, which 
would arise, for instance, from some fundamental failure of the post-Newtonian expan- 
sion to approximate the physics of black holes at high order Simply, they reflect some 
inconsistency, of mathematical origin, in the Hadamard regularization scheme, when it is 
applied to the computation of certain integrals at the 3PN order. Alternatively, one can 
say that this regularization, when "literally" pushed to its maximum (in the way proposed 
in [14,15]), reveals some "incompleteness" in making physical predictions, which can or 
cannot be removed by external physical arguments. Fortunately, we shall see that the ambi- 
guity constant (1) can be resolved once one disposes of the appropriate mathematical tools 
for performing the regularization. 

An improved version of the Hadamard regularization, defined in [14,15], is based on: (i) 
Systematic use of "partie-finie" pseudo-functions to represent the functions in the problem 
which are singular at the location of the particles; (ii) Specific distributional derivatives 
generalizing those of the standard distribution theory [3] in order to differentiate the latter 
pseudo-functions; (iii) "Lorentzian" way of performing the regularization, defined by the 
Hadamard partie finie calculated within the Lorentzian rest frame of the particles. We shall 
refer to that regularization [14,15] as the "extended" Hadamard (EH) one. 

The EH regularization constitutes the first step of a complete calculation of the 3PN 
equations of motion [14,15]. The second step, aimed at removing the incompleteness A, 
consists of going to d-dimensional space and using complex analytic continuation in d, in 
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what is known as the dimensional regularization (henceforth abbreviated as "dimreg")-^ 
For the moment it is not possible to derive the 3PN equations of motion in any d dimen- 
sions, i.e. not necessarily of the form d — 3 + e, where e ^ 0. This is why one still 
has to rely on the 3-dimensional calculation of the equations of motion by means of EH 
regularization. This second step (dimensional continuation in d) has already been achieved 
in the context of the 3PN Hamiltonian in ADM coordinates, with result [17] 

w« = . (2) 

In the present contribution we describe our own application of dimreg (so to say "on the 
top" of Hadamard's regularization) to the derivation of the 3PN equations of motion, in the 
framework of harmonic coordinates, based on recent work in collaboration with Damour 
and Esposito-Farese [18]. 



2. Hadamard regularization of Poisson-like integrals 

Let us start by giving some reminders of the way we compute the Hadamard regularization 
of some potentials having the form of Poisson or Poisson-like integrals. Let F(x) be a 
smooth function on M?, except at the value of two singular points yi and y2, around which 
it admits some Laurent expansions of the type (VAT e N) 

F(x)= ^f/p(ni)+o(rn, 0) 

PQ<p<N 

where n = |x — yi | — > 0, and the i/p(ni)'s denote the coefficients of the various powers 
of ri, which are functions of the positions and velocities of the particles, and of the unit 
direction rii = (x — yi)/ri of approach to singularity 1 (we have also the same expansion 
corresponding to the singularity 2). The powers of r\ are relative integers, p € Z, bounded 
from below by some typically negative po depending on the function F. 

We shall discuss the prescription (taken in [14]) to define the "value at x' = yi" of 
the singular Poisson integral -P(x') of the source function F{'k). The potential i^(x') is 
defined, at any field point x' different from the singularities, in the sense of the Hadamard 
partie-finie (Pf) of an integral, i.e. 

nx') = -i^pf.,./^nx). (4) 

This "partie finie" involves two constants, si and S2, which parametrize some logarithmic 
terms, and are associated with the characteristics of the regularizing volumes around the 
two particles, which have been excised from in order to define the partie finie by means 
of the limit, when the size of these volumes tends to zero, of the integral external to the 
volumes. 



^Dimreg was invented as a mean to preserve the gauge symmetry of perturbative quantum gauge 
field theories [16]. 
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The value at x' = yi of the function P(x') is defined by the Hadamard partie finie 
in the singular limit x' yi, given as usual by the angular average of the coefficient 
of the zeroth power of r'^ = |x' — yi| when — > 0. Notice first that ^'(x') does not 
admit an expansion when — > of the same type as in Eq. (3), since it involves also 
a term proportional to the logarithm of r[ . Thus we shall have, rather than a power- hke 
expansion, 



P(x') 



E 

p'o<V<N 



I P 



gp{n[) + hpin[)\nr[ 



+ o{ri ), 



(5) 



where the coefficients igp and ihp depend on the angles n'^, and also on the constants si 
and S2, in such a way that when combining together the terms in (5) the constant r[ always 
appears in "adimensionaUzed" form hke in \n{r[ /s\). Then we define the Hadamard partie 
finie at point 1 in the standard way (taking the spherical average of the zeroth-order power 
of r'l), except that we include the contribution linked to the logarithm of r'l, which is 
possibly present into that coefficient. More precisely, we define 



(P)i = (5o) + (-^o)lnr'i 



(6) 



where the brackets denote the angular average, over the solid angle element df2(n'j ) on the 
unit sphere. Let us emphasize that in (6) we have introduced in fact a new regularization 
scale denoted r^, which can be seen as some "small" but finite cut-off length scale [so 
that Inr^ in Eq. (6) is a finite, but "large" cut-off dependent contribution]. To compute 
the partie finie one must apply the definition (6) to the Poisson integral (4), which involves 
evaluating correctly the angular integration therein. The result, proved in Theorem 3 of 
[14], reads 



(7) 



The first term is simply the value of the potential at the point 1, namely P{yi), which 
would in fact constitute a "naive" way to implement the regularization, but would not 
yield 3PN equations of motion compatible with basic physical properties such as energy 
conservation. The supplementary term makes the partie finie to differ from the naive guess 
P(yi) in a way which was found to play a significant role in the computations of [10]. The 
apparent dependence of the result (7) on the scale si is illusory. The si -dependence of the 
R.H.S. of Eq. (7) cancels between the first and the second terms, so the result depends only 
on the constants r'l and S2, and we have in fact the following simpler rewriting of (7), 



— nx)-</-2) 



(8) 



Similarly the regularization performed at point 2 will depend on r2 and Si, so that the 
binary's point-particle dynamics depends on four (a priori independent) length scales r^, 
S2 and r'2, si. Because we work at the level of the equations of motion (instead of, say, the 
Lagrangian), many of the terms we shall need are in the form of the gradient of a Poisson 
potential. For the gradient we have a formula analogous to (8) and given by Eq. (5.17a) of 
[14], namely 
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{diP)i = -^Pf«i,<^2 / d^^'^Fi^) + In (^j-^ {n\ f _i) (9) 

= -^Pf^;,«./'^'^^^W' (10) 

where we have taken into account (in the rewriting of the second line) the always correct 
fact that the constant si cancels out and gets "replaced" by r[. 

We must also treat the more general case of potentials in the form of retarded integrals, 
but because we shall have to consider (in Section 3. below) only the difference between 
dimreg and the Hadamard regularization, it will turn out that the first-order retardation 
(IPN relative order) is sufficient for this purpose. Here we are not interested in radiation- 
reaction effects, so we shall use the symmetric (half-retarded plus half-advanced) integral. 
At the IPN order we thus have to evaluate 

i?(x') = P(x') + ^Q(x') + o(^^, (11) 

where -P(x') is given by (4), and where Q(x') denotes (two times) the double or "twice- 
iterated" Poisson integral of the second-time derivative, still endowed with a prescription 
of taking the Hadamard partie finie, namely 

Q(x') = -^Pf.i,.. j c^'x |x - x'|a|i^(x) . (12) 

In the case of Q(x') the results concerning the partie finie at point 1 were given by 
Eqs. (5.16) and (5.17b) of [14], 

= -^Pfri,«. j rf'xria,2F(x) + _4) , (13) 
(5iQ)i = ^Pfr;,.. j d'^n\d^F{^) + i(nl k -3> , (14) 

where the i/cp's denote the analogues of the coefficients i/p, parametrizing the expansion 
of F when r[ 0, but corresponding to the double time-derivative dfF instead of F. 

There is an important point concerning the treatment of the repeated time derivative 
dfF{x) in Eqs. (14). As we are talking here about Hadamard-regularized integrals (which 
excise small balls around both yi and y2), the value of dfF{x) can be simply taken in the 
sense of ordinary functions, i.e., without including eventual "distributional" contributions 
proportional to S{x — yi) or 5(x — y2) and their derivatives. However, we know that such 
terms are necessary for the consistency of the calculation (without them, for instance, the 
calculation would be incorrect already at the 2PN order). In EH regularization, there is a 
specific prescription for the distributional derivative which is issued from the generaUzed 
framework of [14]. In dimreg we shall use simply the standard Schwartz distributional 
derivatives in d dimensions. [As it turns out, the Schwartz derivatives yield some ill- 
defined (formally infinite) expressions in 3 dimensions — this is why a generalization of 
the Schwartz distributional derivative defined in [14] was found to be necessary — but 
the latter expressions are proved to be rigorously zero when computed in d dimensions.] 
Therefore, in our computation of the difference between dimreg and Hadamard regulariza- 
tion (next Section), we must also include the difference between the different prescriptions 
for the distributional derivatives in d and in 3 dimensions. We refer to [18] for the details. 
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3. Difference between the dimensional and Hadamard regularizations 

In dimreg the computation of the regularized value of Poisson or Poisson-like integrals is 
very simple [18]. First of all, the generalization of the function F \n d dimensions will be 
some F'''' which admits when ri ^ a more complicated expansion, namely (ViV € N) 



F('')(x)= Y. ^r'V^:^(ni)+o(rf), (15) 

po<p<N ^ 

go<q<qi 



where p and q are relative integers (p, q G Z), whose values are limited by some po, 
qo and qi as indicated. The expansion (15) differs from the corresponding expansion in 3 
dimensions, as given in Eq. (3), by the appearance of integer powers of rf where we denote 
e = d — 3. The coefficients i/pf, depend on the unit vector ni in d dimensions, on the 
positions and coordinate velocities of the particles, and also on the characteristic length 
scale £o of dimensional regularization. The latter can be introduced into the formalism by 
saying that the constant G in the R.H.S. of the Einstein field equations is related to Gn, the 
usual Newton constant in 3 dimensions, by G = GnIq- Because F^"*' F when d ^ 3 
we necessarily have some constraint on the coefficients i/pf, so that we are in agreement 
with the expansion (3) in this limit. 

Consider now the Poisson integral of F^'^\ in d dimensions, given by the standard 
Green's function for the Laplace operator in d dimensions, namely 

P«^)(x') = A- [F('^)(x)] . - A / r^^^i^('')(x) , (16) 



where k is related to the usual Eulerian F-function by 



3 



k=-W-- (17) 



TT 2 



To evaluate the Poisson integral at the singular point x' = yi is quite easy to do in dimreg, 
because the nice properties of analytic continuation allow simply to get [F'('')(x')]x'=yi 
by replacing x' by yi into the explicit integral form (16). So, we simply have, 

p^''\y^) = -ll^Fi^\^). (18) 

Similarly, for the twice iterated Poisson integral, and the relevant gradients of potentials. 



^We have limd_,3 k = 1. Notice the following connection to the volume of the sphere with d — 1 
dimensions {i.e., embedded into Euclidean d-dimensional space): 

47r 

kfld-i = -, — - ■ 
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Q^'Hyi) = -4^(^rrf) / rf'xr^'^a,2i^('^)(x) , (19) 

diP^'Kyi) = -^^^ / '^'^x ^F^''\^) , (20) 

diQ^^Hyi) = ^J d''xn\rl-''d^F^''\K) . (21) 



The main technical step of our strategy will then consist of computing the difference 
between the d-dimensional Poisson-type potentials (18)-(21), and their "pure Hadamard- 
Schwartz" 3-dimensional counterparts, given by expressions such as (8). By pure 
Hadamard-Schwartz (pHS) we mean in some sense the "core" of the Hadamard regular- 
ization, i.e. merely based on the usual notion of the partie finie of a singular function or a 
divergent integral, but without the improvements brought about by the EH regularization 
(see [18] for more details). For instance, the computations of Section 2. above belong to 
the pHS regularization, but the special treatment of distributional derivatives in three di- 
mensions is specific to the EH regularization. Given the results (P)i and P^'^''(yi) of the 
two regularizations, denoting the difference by means of the script letter V, we thus pose 

1?P(l) = p('^)(yi)-(P)i. (22) 

That is, DP(1) is what we shall have to add to the pHS result in order to get the correct d- 

dimensional result. Note that, in this paper, we shall only compute the first two terms, 
a_i£^^ + oq + 0{e), of the Laurent expansion of T>P{1) when e ^ 0. [We leave 
to future work an eventual computation of the rf-dimensional equations of motion as an 
exact function of the complex number d.\ This is the information we shall need to fix 
the value of the parameter A. As we shall see, the difference 'DP(l) comes exclusively 
from the contribution of poles oc 1/e (and their associated finite part) in the rf-dimensional 
calculation. Here we simply state the result without proof (see [18] for details). We obtain 
the following closed-form expression for the difference [valid up to the neglect of higher- 
order terms 0(e)], 

^^<" — ^ G+4>»K-ll)(/L-L> (23, 

which constitutes the basis of all the practical calculations in the work [18].* Here we still 
use the bracket notation to denote the angular average, but now performed in d dimensions, 



''with the same notation as in [18] the multipole expansion in d dimensions reads as 

2-d (~)^ a ( 1 A I L 
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i.e. over the solid-angle element dVLd-i associated with the {d — l)-dimensional sphere. 
Notice that (23) depends on the two "constants" Inr'j^ and lns2. As we shall check, these 
In r'^ and In S2 will exactly cancel out the same constants present in the pHS calculation, so 
that the dimreg acceleration will be finally free of the constants and S2. Note also that 

the coefficients ifp], and in d dimensions depend on the length scale £0 associated 
with dimreg. Taking this dependence into account one can verify that and 82 in (23) 
appear only in the combinations ln(r5^/£o) and ln(s2/£o)- 

Let us give also the formula for the difference between the gradients of potentials, i.e. 

Vd,P{l) = d^P^''\yi)-{d,P)i, (24) 

which is readily obtained by the same method. We have 



VdiP{l) = - i ^ [i + elnrl) (nl / (25) 

Formulae (23) and (25) correspond to the difference of Poisson integrals. But we have 
already discussed that we need also the difference of inverse d' Alembertian integrals at the 
IPN order. To express as simply as possible the IPN-accurate generalizations of Eqs. (23) 
and (25), let us define two junctionals H and Hi which are such that their actions on any 
(/-dimensional function f ("^^ is given by the R.H.S.'s of Eqs. (23) and (25), i.e., so that 

VP{1) =n[F'-'^^ , (26) 

pa,;P(i) = [f'^'^^ . (27) 

The difference of IPN-retarded potentials and gradients of potentials is denoted 

VRil) = R(^\yi)-{R)i, (28) 

Vd,R{l) = d,R<^''\yi) - {d,R)i , (29) 

where in 3 dimensions the potential R{x') is defined by Eq. (11) and the regularized values 
(i?)i and {diR)i follow from (8), (10), (14), and where in d dimensions R^'^''{yi) and 
5ii?*^'^)(yi) are consequences ofEqs. (18)-(21). With this notation we now have our result, 
that the difference in the case of such IPN-expanded potentials reads in terms of the above 
defined functionals H and Hi as 



VR{1) = n 



VdiR{l) = Hi 



^2c^{A-dy'^ 
2c2(d-2) * 



-i,{k-.) + o[^), (30) 

-i^(-1^3>^-0(^)• (31) 



These formulae involve some "effective" functions which are to be inserted into the func- 
tional brackets of H and Hi. Beware of the fact that the effective functions are not the same 
in the cases of a potential and the gradient of that potential. Note the presence, besides the 
main terms H[- ■ ■] and Hi[- • •], of some extra terms, purely of order IPN, in Eqs. (30)- 
(31). These terms are made of the average of some coefficients ikp of the powers in 
the expansion when ri ^ of the second-time-derivative of F, namely dfF. They do not 
seem to admit a simple interpretation. They are important to get the final correct result. 
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4. Dimensional regularization of the equations of motion 

We outline next the way we obtain from the previous computation of the "difference" the 
3PN equations of motion in dimreg, and show how they are physically equivalent to the 
EH-regularized equations of motion. We start from the end results of [10] for the 3PN 
acceleration of the first particle, say af^, depending on the two arbitrary length scales 
r'l and r'2 (appearing when regularizing Poisson-like integrals in Section 2.), and on the 
"ambiguity" parameter A. Explicitly, we define 



Here the acceleration is considered as a function of the two masses mi and 1712, the relative 
distance yi — y2 = ?'i2ni2 (where ni2 is the unit vector directed from particle 2 to particle 
1), the two coordinate velocities Vi and V2, and also, as emphasized in (32), the parameter 
A as well as two regularization length scales r[ and r'2. The latter length scales enter the 
equations of motion at the 3PN level through the logarithms ln(ri2/?'^) and In(ri2/r2). 
They come from the regularization as the field point x' tends to yi or y2 of Poisson-type 
integrals (see Section 2. above). The length scales r'^, r'2 are "pure gauge" in the sense that 
they can be removed by the effect induced on the world-lines of a coordinate transformation 
of the bulk metric [10]. 

On the other hand, the dimensionless parameter A entering the final result (32) corre- 
sponds to genuine physical effects. It was introduced by requiring that the 3PN equations 
of motion admit a conserved energy (and more generally be derivable from a Lagrangian). 
This extra requirement imposed two relations between the two length scales r'^, r'2 and the 
two other length scales si, S2 entering originally into the formalism, namely the constants 
si and S2 parametrizing the Hadamard partie finie of a Poisson integral as given by Eq. (4) 
above. Recall that s\ and S2 are associated with the characteristics of the two regularizing 
volumes (notably their shape) around the singularities, which are excised in order to define 
the Hadamard partie finie of a divergent integral. The latter relations were found to be of 
the form 



(and 1 ^ 2), where the so introduced single dimensionless parameter A has been proved 
to be a purely numerical coefficient {i.e. independent of the two masses). It is often con- 
venient to insert Eq. (33) into (32) and to reexpress the acceleration of particle 1 in terms 
of the original regularization length scales entering the Hadamard regularization of ai, 
which were in fact r'^ and S2 [as shown, for instance, in Eq. (8)]. Thus we can consider 
alternatively 



where the regularization constants are subject to the constraints (33) [we can check that 
the A-dependence on the R.H.S. of (34) disappears when using Eq. (33) to replace r'2 as a 
function of S2 and A]. 

The strategy followed in [18] consists of two steps. first step consists of subtracting 
all the extra contributions to Eq. (32), or equivalently Eq. (34), which were specific conse- 
quences of the EH regularization defined in [14, 15]. As has been detailed in [18], there are 



a?"" [A; r'l, r^] = R.H.S. of Eq. (7.16) in Ref. [10] . 



(32) 




(33) 



af[r;,,S2]^af[A;r'i,r^(s2,A)], 



(34) 
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seven such extra contributions 5^sli, A = \. ■ ■ ■ ,7 . Subtracting these contributions boils 
down to estimating the value of ai that would be obtained by using a "pure" Hadamard 
regularization, together with Schwartz distributional derivatives, which is what we have 
already called the "pure Hadamard-Schwartz" (pHS) regularization. Such a pHS accelera- 
tion was in fact essentially the result of the first stage of the calculation of ai, as reported 
in the (unpublished) thesis [19]. It is given by 

ar[ri,S2] =af[ri,S2]-E'^^^i- (35) 

A=l 

The second step of our method consists of evaluating the Laurent expansion, in powers of 
e = d — 3, of the difference between the dimreg and pHS (3-dimensional) computations 
of the acceleration ai. As we have seen in Section 3. this difference makes a contribution 
only when a term generates a pole ^ 1/e, in which case dimreg adds an extra contribution, 
made of the pole and the finite part associated with the pole [we consistently neglect all 
terms 0{e)]. One must then be especially wary of combinations of terms whose pole 
parts finally cancel ("cancelled poles") but whose dimensionaUy regularized finite parts 
generally do not, and must be evaluated with care. We denote the above defined difference 

Pai = Vai [e, Iq; r[, S2] = Va^ [s, Iq; A; r'l, r^] . (36) 

It depends both on the Hadamard regularization scales r[ and S2 (or equivalently on A and 
r[, r'2) and on the regularizing parameters of dimreg, namely e and the characteristic length 
Iq. It is made of the sum of all the individual contributions of the Poisson or Poisson-like 
integrals as computedin Section 3. above [e.g. Eqs. (23)and(25)]. Finally, our main result 
will be the explicit computation of the e-expansion of the dimreg acceleration as 

af"^^[£,4] = af\[,S2] + 1?ai[£,4; r'l, S2] . (37) 

With this result in hands, we have proved [18] two theorems. 

Theorem 1 The pole part oc 1/e of the dimreg acceleration (37), as well as of the metric 
field gfiu{x) outside the particles, can be re-absorbed (i.e., renormalized away) into some 
shifts of the two "bare " world-lines: ya — > ya + ^a, with, say, oc 1/e ("minimal 
subtraction"; MS), so that the result, expressed in terms of the "dressed" quantities, is 
finite when e ^ 0. 

The situation in harmonic coordinates is to be contrasted with the calculation in ADM-type 
coordinates within the Hamiltonian formalism [17], where it was shown that all pole parts 
directly cancel out in the total 3PN Hamiltonian (no shifts of the world-lines were needed). 
The central result of the paper is then as follows. 

Theorem 2 The "renormalized" (finite) dimreg acceleration is physically equivalent to 
the EH-regularized acceleration (end result ofRef. [10]), in the sense that there exist some 
shift vectors ^i(s,£o;r[) and ^2(£, ^o! ?'2)> suchthat 

^i^lK r[,r',] = lim [af^'^'^^[s, io] + 5^(e,^o;r-l,ri) ai] (38) 
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(where 5^ ai denotes the effect of the shifts on the acceleration if and only if the hereto- 
fore unknown parameter A entering the harmonic-coordinates equations of motion takes 
the value 



The precise shifts ^„ (e) needed in Theorem 2 involve not only a pole contribution oc 1/e, 
which defines the "minimal" (MS) shifts considered in Theorem 1, but also a finite contri- 
bution when e — > 0. Their explicit expressions read: 



^1 = 



11 



2 

mi 




327 
1540 



ajvi and 1 2 , 



(40) 



where Gat is the usual Newton's constant, ajvi denotes the acceleration of the particle 1 
(in d dimensions) at the Newtonian level, and q = Ane'-" depends on the Euler constant 
G = 0.577- ••. 

An alternative way to phrase the result (38)-(39), is to combine Eqs. (35) and (37) in 
order to arrive at 



lim 



ai 



(41) 



Under this form one sees that the sum of the additional terms (5^ai differs by a mere 
shift, when and only when X takes the value (39), from the specific contribution Vai, 
which comes directly from dimreg. Therefore one can say that, when A = — ggfg, the 
EH regularization [14,15] is in fact (physically) equivalent to dimreg. However tiie EH 
regularization is incomplete, both because it is a priori unable to determine A, and also 
because it necessitates some "external" requirements such as the imposition of the link 
(33) in order to ensure the existence of a conserved energy — and in fact of the ten first 
integrals linked to the Poincare group. By contrast dimreg succeeds automatically (without 
extra inputs) in guaranteeing the existence of the ten conserved integrals of the Poincare 
group, as already found in Ref. [17]. 

In view of the necessary link (1) provided by the equivalence between the ADM- 
Hamiltonian and the harmonic-coordinates equations of motion, our result (39) is in perfect 
agreement with the result = obtained in [17]. [One may wonder why the value of A 
is a compUcated rational fraction while Wg is so simple. This is because u)s was introduced 
precisely to measure the amount of ambiguities of certain integrals, while, by contrast, A 
has been introduced as the only possible unknown constant in the link between the four 
arbitrary scales r'i,r2, Si,S2 (which has a priori nothing to do with ambiguities of inte- 
grals), in a framework where the use of the EH regularization makes in fact the calculation 



^When working at the level of the equations of motion (not considering the metric outside the 
world-lines), the effect of shifts can be seen as being induced by a coordinate transformation of the 
bulk metric as in Ref. [10]. 



Pramana - J. Phys., Vol. 53, No. 6, December 1999 



11 



Luc Blanchet 



to be unambiguous.] Besides the confirmation of the value of ojs or A, this result pro- 
vides a confirmation of the consistency of dimreg, because our explicit calculations are 
entirely different from the ones of [17]: We use harmonic coordinates (instead of ADM- 
type ones), we work at the level of the equations of motion (instead of the Hamiltonian), we 
use a different form of Einstein's field equations and we solve them by a different iteration 
scheme. Our result is also in agreement with the recent finding of Refs. [20,21] (see also 
[22]), where the 3PN equations of motion are derived in harmonic gauge using a "surface- 
integral" approach, aimed at describing extended relativistic compact binary systems in the 
strong-field point particle limit. 



5. Equations of motion of circular-orbit compact binaries 



From a practical point of view, the determination of the value of A allows one to use the 
full 3PN accuracy in the analytical computation of the dynamics of the last orbits of binary 
systems [23,24]. We assume a circular orbit since most inspiralling compact binaries will 
have been circularized at the time when they enter the frequency bandwidth of the detectors 
LIGO and VIRGO. In the case of circular orbits — apart from the gradual 2.5PN radiation- 
reaction inspiral — the quite complicated equations of motion, Eq. (7.16) in Ref. [10], 
simplify drastically. 

We translate the origin of coordinates to the binary's center-of-mass by imposing that 
the binary's center-of-mass vector, deduced from the Lagrangian formulation of the 3PN 
equations of motion, is zero (see e.g. Ref. [25]). Then, in the center-of-mass frame, the 
relative acceleration ai2 = ai — a2 of two bodies moving on a circular orbit at the 3PN 
order is given by 



ai2 



-w^yi2 



32 



-Vl2 



12 



(42) 



where yi2 = yi — y2 is the relative separation (in harmonic coordinates) and lo denotes 
the angular frequency of the circular motion; the second term in Eq. (42), opposite to the 
velocity vi2 = Vi — V2, is the 2.5PN radiation reaction force. In (42) we have introduced, 
in addition to the total mass m = mi+ m2, the symmetric mass ratio 



m\m2 



(43) 



which is generally very useful because of its interesting range of variation Q < v < \, 
with 1/ = i in the case of equal masses, and i/ — » in the "test-mass" limit for one of the 
bodies. The main content of the 3PN equations (42) is the relation between the frequency 
uj and the orbital separation ri2, that we find to be given by the 3PN-generaUzed "Kepler" 
third law [9,10] 



''SPN 



12 



41 



(44) 



+ -10- 



75707 41 



840 



22ln 
64 \r'a 



19 
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in which we employ, in order to display the successive post-Newtonian corrections, the 
post-Newtonian parameter [of the order of 0(l/c^)] 



7: 



Gm 



(45) 



The acceleration (42)-(45) is entirely specified at the 3PN order except for some unphysical 
gauge freedom, parametrized by the length scale Tq appearing in Eq. (44), which is nothing 
but the "logarithmic barycenter" of the two gauge-constants and entering the end 
results of [10], i.e. 



mi TO2 , 

m Tf, = in Ti H in rn 

m m 



(46) 



As for the binary's energy (in the center-of-mass frame), it is readily obtained from the 
circular-orbit reduction of the conserved energy associated with the 3PN Lagrangian in 
harmonic coordinates [12]. We find 



i?3PN 



1 



235 



7 1 

-4 + 4^)7 + 
46031 123 



7 49 



22 



2240 



64 



In 



ri2 



1 



7 

27 
— I 
32 



(47) 



64 



This expression is that of a physical observable E, however it depends on the choice of a 
coordinate system, because it involves the post-Newtonian parameter 7 defined from the 
harmonic-coordinate separation ri2. But the numerical value of E should not depend on 
the choice of a coordinate system, so E must admit a frame-invariant expression, the same 
in all coordinate systems. To find it we re-express E with the help of a frequency-related 
parameter x instead of the separation-related parameter 7 [this is always a good thing to 
do] . We define a; to be, like for 7, of the order of O ( 1 /c^ ) by posing 



/ Gmui 



2/3 



(48) 



Then we readily obtain the expression of 7 in terms of x at 3PN order, 

65 



73PN 



= .{l+( 



+ 1 



V 

L 

3 

2203 
2520 



x+ ( 1 
41 



12 



V X 



2 22 

TT m 

192 3 



r\i \ 
r'o) 



(49) 



229 o 1 \ ^ 

u -\ H x'^ 

36 81 



that we substitute back into Eq. (47), making all appropriate post-Newtonian re- 
expansions. As a result we gladly discover that the logarithms together with their asso- 
ciated gauge constant have cancelled out. Therefore our final result is 



-E'SPN — 



, ,3 1 \ / 27 19 1 2\ 2 



675 



34445 205 



96 



576 



155 
-96' 



35 
5184' 



x^ 
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In the test-mass limit ^ 0, we recover the energy of a particle with mass /i = m.v in 
a Schwarzschild background of mass m, i.e. -Etest = /^c^ [(1 ^ 2x)(l — 3a;)^^/^ — l], 
when developed to 3PN order. Of course, the subtleties we have discussed, linked with the 
self-field regularization, disappear in the test-mass Umit, but, interestingly enough, they 
affect only the term proportional to u in the 3PN coefficient of Eq. (50); the terms propor- 
tional to and in Eq. (50) have been found to be "complete" in EH regularization. 



NUMERICAL 
* 

3PN 
* 2PN 



0.1 0.2 0.3 0.4 0.5 0.6 



Figure 1. Results for Eico in terms of wico for equal-mass binaries (v = i). The 
asterisk marks the result calculated by numerical relativity. 

Finally let us compute the innermost circular orbit (ICO) of point-particle binaries 
through 3PN order, following [24]. The ICO is defined as the minimum, when it exists, of 
the binary's circular-orbit energy function (50). In particular, we do not define the ICO as a 
point of dynamical (general-relativistic) unstability. [See Section 6 of [25] for a discussion 
of the dynamical unstability in the post-Newtonian framework.] In Fig. 1 we plot Eiqq 
versus wico in the case of equal masses {v = \), and compare the values with the recent 
finding of numerical relativity, obtained by means of a sequence of quasi-equiUbrium con- 
figurations under the assumptions of helical symmetry and conformal flatness [26,27]. As 
we can see the 2PN and 3PN points are rather close to each other and to the numerical 
value. However, the IPN approximation is clearly not precise enough, but this is not very 
surprising in the highly relativistic regime of the ICO where the orbital velocity reaches 
v/c ~ {GmujiGo/ c^Y^^ ~ 0.5. A striking fact from Fig. 1 is that the post-Newtonian 
series seems to "converge well", but actually the series could be only asymptotic (hence 
divergent), and, of course, still give excellent results provided that the series is truncated 
near some optimal order of approximation. 

Our conclusions, therefore, are that (1) the post-Newtonian approximation is likely to be 
vaUd and quite accurate in the regime of the ICO (in the equal-mass case), and (2) it is in 
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good agreement with the result of numerical relativity. Note that the conclusion (1) contra- 
dicts some earlier prejudices about the slow convergence of the post-Newtonian approxi- 
mation (see e.g. Ref. [4]). Furthermore, our computations are based on the standard post- 
Newtonian expansion, without using any resummation techniques such as Pade approxi- 
mants and/or effective-one-body method. For recent comparisons of the post-Newtonian 
and numerical calculations in the regime of the ICO, including finite-size effects appropri- 
ate to neutron-star binaries, see [28,29]. 
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